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a^ ■ ABSTRACT 

o : 

O . The luminosity function of white dwarfs is a powerful tool for studies of the evolution 

O I and formation of the Milky Way. The (theoretical) white dwarf cooling sequence provides 

,-^ ■ a useful indicator of the evolutionary time scales involved in the chronometry and star 

formation history of the galactic disk, therefore, intrinsically faint (& old) white dwarfs 
O I in the immediate solar neighborhood can be used to determine an upper limit for the 

^ I age of the galactic disk. 

cd ■ Most determinations of the white dwarf luminosity function have relied on the use 

of Schmidt's 1/Knax (Schmidt 1975) method for magnitude and proper-motion selected 
k>< I samples, the behavior of which has been demonstrated to follow a minimum variance 

^ I maximum-likelihood pattern for large samples. Additionally, recent numerical simula- 

tions have also demonstrated that the 1/Knax provides a reliable estimator of the true 
LF, even in the case of small samples (Wood and Oswalt 1998, and Garcia-Berro et 
al. 1999). However, the conclusions from all these previous studies have been based on 
noise-free data, where errors in the derived LF have been either assumed to follow a 
Poisson distribution (valid only for large samples), or where other simple estimates for 
the uncertainties have been used. 

In this paper we examine the faint-end {My > +14) behavior of the disk white 
dwarf luminosity function using the 1/Knax method, but fully including the effects of 
realistic observational errors in the derived luminosity function. We employ a Monte 
Carlo approach to produce many different realizations of the luminosity function from a 
given data set with pre-specified and reasonable errors in apparent magnitude, proper- 
motions, parallaxes and bolometric corrections. These realizations allow us to compute 
both a mean and an expected range in the luminosity function that is compatible with 
the observational errors. 



We find that current state-of-the art observational errors, mostly in the bolometric 
corrections and trigonometric parallaxes, play a major role in obliterating (real or ar- 
tificial) small scale fiuctuations in the luminosity function. We also find that a better 
estimator of the true luminosity function seems to be the median over simulations, rather 
than the mean. When using the latter, an age of 10 Gyr or older can not be ruled out 
from the sample of Leggett, Ruiz, and Bergeron (1998). 

Subject headings: Galaxy: formation — solar neighborhood — stars: white dwarfs — 
methods: data analysis — methods: numerical — surveys — catalogs 



Introduction 



2. The luminosity function 



The luminosity function (LF) for (disk) white dwarfs 
(WDs) is an important observational tool to guide our 
understanding of the evolution and fate of intermediate- 
mass systems, as well as for the more general Galactic 
structure problems of the age of the Galactic disk, as 
first pointed out by D'Antona and Mazzitelli (1978), 
and the star-formation history of the solar neighbor- 
hood (Isern et al. 1995). Indeed, WDs evolutionary 
time-scales represent a useful tool for constraining the 
age of the disk of our Galaxy. The existence of an 
abrupt fall-off in the observed WD LF (see Liebert, 
Dahn and Monet (1988), Leggett, Ruiz and Berg- 
eron (1998)) has been interpreted as an indication 
of the finite age of the Galactic disk (Winget et al. 
1987). By fitting the observations with theoretical 
WD LFs, this interpretation has been quantitatively 
explored by various investigators (Wood 1992, and 
references therein). Not considered in these studies 
has been an analysis of the effects of observational 
errors on the observationally derived LF. Recently, 
Wood and Oswalt (1998) and Garcia-Berro et al. 
(1999), have performed a very comprehensive set of 
numerical simulations studying the effects of different 
star-formation rates, IMF, and kinematical prescrip- 
tions for (the progenitors of) WDs on the derived the- 
oretical WD LF. However, all their simulations were 
performed under the assumption that observational 
quantities were noise-free. 

In this paper we attempt to realistically quantify 
and characterize the effects of observational errors on 
the derived LF for WDs by comparing a LF derived 
assuming no errors at all, with LFs using errors on 
the various basic observational parameters, applying 
the 1/Vmax method. 

In Section 2 we present a general discussion of the 
importance and scope of determining the LF for WDs. 
In Section 3 a brief description of the l/Knax method 
is given, while Section 4 describes recent determina- 
tions of the WD LF using this method. In Section 5 
we describe the basis of our numerical simulations. 
Sections 6 and 7 present the results of our simulations 
in terms of the global uncertainties in the WD LF us- 
ing current data and the effects of individual sources 
of errors, respectively. Section 8 outlines our conclu- 
sions. 



The classical method for determining the LF of 
magnitude and proper-motion selected samples is that 
proposed by Schmidt (1975). This method, called 
the 1/Vmaxj stems from a generalization of a method 
proposed earlier by Schmidt (1968) for magnitude- 
limited samples. The method assumes that the LF 
does not change (evolves) as a function of distance 
from the observer, and that the sample is homoge- 
neously distributed in space. The 1/Vmax method 
computes the LF by weighting the contribution of 
each observed point by the equivalent volume where 
that particular object could have been observed under 
the pre-specified survey constraints. Felten (1976) has 
shown that the 1/Vinax method for magnitude-limited 
samples is a minimum variance maximum-likelihood 
estimator and that, for small absolute magnitude bins 
and very large samples (> 200 objects per bin, not 
the case of WDs), it provides a reliable way of esti- 
mating the parent (true) LF. Several modifications 
have been proposed to the method in the case of 
magnitude-limited samples (Davis and Huchra 1982, 
Bales 1993), most notably one that allows the combi- 
nation of different samples coherently (Avni and Bah- 
call 1980). However, the basic scheme to determine 
the LF of magnitude and proper- motion selected sam- 
ples has remained unchanged, and few and limited nu- 
merical simulations have been carried-out to explore 
the robustness and possible biases that the original 
method might have when dealing with complete, but 
small, kinematically selected samples (Wood and Os- 
wah 1998, Garcia-Berro et al. 1999). 

Because the spatial density of WDs is rather small 



(about 3.4 X 10 -^ stars/pc"^ down to Mv 



-16.75), 



it is important to ensure that the method used to de- 
termine its LF is either free from biases, or that they 
can be at least reliably corrected. Also, it is impor- 
tant to understand the effects of the kinematic selec- 
tion on the resulting LF. For this purpose. Wood and 
Oswah (1998) and Garcia-Berro et al. (1999), have 
performed extensive numerical simulations by creat- 
ing fake catalogues of WDs in the solar neighborhood 
from a pre-specified LF and a kinematical description. 
Their (predicted) mock catalogues are an approxima- 
tion to true catalogues with similar selection biases in 
apparent magnitude and total proper-motion. These 
mock catalogues are then passed to the 1/Vmax from 
which a LF is predicted. This LF is then compared 
to a range of input LFs, with different catalogue con- 



straints, and selection effects. The main results from 
the Wood and Oswalt (1998) work are 1) the l/Vmax 
method provides robust estimates of the true local 
space density, 2) the age of the galactic disk must be 
considered uncertain by about 15% for the currently 
available sample sizes, and 3) the bright-end of the de- 
rived LF shows substantial deviations from the input 
functions, suggesting that is is difficult to derive vari- 
ations in the recent star formation history of the disk 
from magnitude and proper- motion selected samples. 
Similarly, Garcia-Berro et al. (1999) also find that 
1) the simulated and observed LFs are in excellent 
agreement, 2) the effect of a scale-height law are im- 
portant, specially at large intrinsic luminosities (i.e., 
the bright end of the LF, which we do not consider 
in this paper), 3) observational errors in the LF are 
well represented by Poisson errors for samples of 200 
stars or more, and 4) the statistical uncertainty in the 
age of the disk is about 1 Gyr, in agreement with the 
findings from Wood and Oswalt. 

Unfortunately, the Wood and Oswalt's and Garcia- 
Berro's et al. simulations have not included the effect 
of realistic observational random errors in the key ob- 
servational quantities, and thus the effect of these er- 
rors on the resulting LF has not been evaluated. This 
is precisely the motivation and scope of this work. Of 
course, of particular interest, is the behavior of these 
simulations with respect to the slope of the faint end 
of the WD LF, which, as mentioned in Section (Q), 
can be used as a constraint for the age of the Galac- 
tic disk. Another point of interest is the level at which 
the detailed shape ( "wiggles" ) on the WD LF are real, 
for a given sample-size, and can be interpreted as a 
consequence of the evolution of WDs as a function of 
cooling age (Diaz-Pinto et al. 1994, Hernanz et al. 
1994). 

3. The 1/Vmax method 

The method proposed by Schmidt (1968, 1975) al- 
lows for a derivation of the LF for a complete and 
spatially uniform sample of stars for which we know 
their apparent magnitudes, parallaxes and (if used in 
the sample selection), proper- motions. We also need 
to know the sample selection (or survey) limits. 

If we have a sample with a lower proper-motion 
limit fj,i and a faint apparent magnitude limit m/, 
the maximum distance rmax over which any star can 
contribute to the sample is given by: 



rmax = mzn[p-i(/i/^0;P"'lO"-'^'"""^^] (1) 

where p is the parallax, fi is the proper-motion, 
and m the apparent magnitude. 

Similarly, if the sample is only complete to an up- 
per proper-motion limit fiu and a bright apparent 
magnitude rub, the minimum distance for inclusion 
would be: 



r^in = maxb-i(/i//i„);p-ilOO-2("'-")] (2) 

Finally, if the sample only covers a fraction f3 of 
the sky, then the maximum volume in which a star 
can contribute to the sample is: 



Knax = ^7r/3(r3 



(3) 



The contribution to the LF from each star in the 
sample is then 1/Knaxj and the LF is calculated by 
adding the l/Vmax values over discrete magnitude in- 
tervals. For more details of the method, the reader is 
referred to Schmidt (1968, 1975). 

As can be seen from the above equations, the con- 
tribution from every star in the sample to the overall 
LF is highly non-linear in terms of the basic observa- 
tional quantities, hence preventing an analytic treat- 
ment of errors, specially in the case of small samples. 
This is even more relevant if we consider that every 
observational point has its own error budget, which 
highlights the need for doing full numerical simula- 
tions of the effect of errors on the derived LF. Indeed, 
Wood and Oswalt (1998) have pointed out that the 
noise properties of the 1/Knax are not well understood 
in the limit of small samples (see also Felten 1976). 

4. Previous WD LF determinations using the 



1/14, 



method 



Liebert, Dahn & Monet (1988, LDM88 thereafter) 
have presented trigonometric parallaxes, optical col- 
ors, and spectrophotometric data for a complete sam- 
ple of intrinsically faint WDs (derived from the Luyten 
Half Second catalogue, LHS, Luyten 1979) in the 
context of a program to determine the faint end 
of the WD LF. Using the classical 1/T4nax method, 
they derived a LF which indicated a downturn near 
log{L/LQ) ~ —4.4, a stellar density of 3 x 10~^ stars 
pc~"^, and a derived age for the disk in the range 



7—10 Gyr. More recently, Legget, Ruiz and Bergeron 
(1998, LRB98 thereafter) gathered new optical and 
near-IR data for the same cool WDs in the LDM88 
sample. Using stellar parameters derived from these 
data and the more refined model atmospheres by 
Bergeron et al. (1995), they re-derived the faint- 
end of the WD LF, also using the 1/T4iax method. 
Comparing their LF with the (then) most recent cool- 
ing sequences by Wood (1995), they derived a rather 
young age for the disk of 8 ± 1.5 Gyr. In both cases, 
the uncertainty on the LF was computed using the 
classical approach of assuming Poisson noise in the 
counts of every absolute (or bolometric) magnitude 
bin, without consideration of the actual observational 
errors for the quantities involved in the LF determina- 
tion. Therefore, only sampling errors were considered, 
but not observational errors. In LRB98, furthermore, 
the authors endeavored to estimate the uncertainties 
in the WD LF by calculating how many stars could 
be thrown into or out of each magnitude bin due to 
errors in the bolometric correction - this is in some 
sense a "precursor" of our numerical simulations (see 
Sect. |.) 

Figure |l| shows a comparison of the LF from LDM88 
and LRB98 as a function of absolute magnitude My, 
adopting a bin size of 0.5 mag similar to that used by 
LDM88 and LRB98. Error bars are Poisson bars, as 
adopted by these authors. We also indicate the LF as 
a function of luminosity, using the same bolometric 
corrections (BCs hereafter) adopted by the authors. 
In both cases (and in what follows of our analysis) 
we have adopted a bright and faint apparent mag- 
nitude limits of Vb — +1, Vf — -1-19, a lower and 
upper proper-motion limits of /i; = 0.8 arcsec yr~^, 
Hu — 10.0 arcsec yr~^, and a fraction of the sky cov- 
ered of /? = 0.5368 (this last value is derived from 
the fact that the LHS catalogue, on which the sample 
is based, covers only the sky north oi d = —20°, and 
avoids the Galactic plane). In the absolute magnitude 
range sampled by these WDs the global normalization 
in the range -1-12.75 < My < -f 16.75 is very similar, 
and equal to p, = 2.46 x 10"^ stars pc"^ for LDM88 
and p^ = 2.54 x 10"^ stars pc^^ ^j. LRB98. 

As we shall see, one of the key ingredients in deter- 
mining a WD LF that could be compared with the- 
oretically derived LFs is the bolometric correction. 
To derive their observational WD LF, LDM88 used 
two extreme bolometric corrections, namely, no cor- 
rection at all, and another one based on the rather 
uncertain model atmospheres available at that time. 



LRB98 on the other hand used not only the latest 
model atmospheres, but they also fitted the detailed 
shape of the theoretical spectrum to the observed op- 
tical and near-IR broad-band colors for every single 
star in the sample, separately (details of the fitting 
technique are given in Bergeron et al. 1997). In this 
last case, errors in effective temperature were derived 
from uncertainties in the fit, while the errors of the 
radius (surface gravity), where derived by propagat- 
ing the uncertainty in the trigonometric parallaxes. 
The bolometric magnitude was then computed using 
Mboi = -2.75 log(L/Lo) + 4.75, with L = AnR^aT^g. 

5. Numerical simulations 

In this section we present the results from our nu- 
merical simulations, fully including observational er- 
rors in all relevant quantities, namely: apparent mag- 
nitude, bolometric corrections, proper-motion, and 
parallaxes. It is assumed that quoted observational 
errors represent the parent standard deviation, and 
that the true value follows a Gaussian distribution 
function with the same standard deviation and mean 
value as that given by the published data. This is 
probably an idealization, but it should provide a bet- 
ter representation of the data than just neglecting the 
observational errors, as it has been so far the case. In 
every single realization of the LF (from now on simply 
called a "simulation" , and usually identified by a se- 
quential integer, or "ID" number), values with a mean 
and dispersion from tabular input quantities taken ei- 
ther from LDM88 or LRB98 are randomly drawn from 
a Gaussian distribution function for all the observ- 
ables (see, e.g., for the case of the simulated parallax, 
Eq. H). The process of generating randomly Gaussian 
distributed values for all the observables of a given 
data point is repeated for all the data points in the 
input sample, thus creating a single simulated sample. 
This simulated sample is then used to construct the 
LF for that particular simulation using the l/Vmax 
method, and the whole process is repeated again to 
create a different simulated sample and its respec- 
tive LF. In the end, various collective values, averaged 
over the simulations, are then produced. In this way, 
it is possible to derive mean, median and quartiles 
for the LF over a given set of simulations, as well as 
other statistical indicators. For example, if <I>*(M) is 
the luminosity function at absolute magnitude M for 
simulation "i" , then the mean-over-simulations LF is 



w=i 



simply given by < $(Af) >= ^^!!p' $*(M)/A^,i 



simul; 



where A^simui is a (pre-specified) number of simulated 



LFs that have been generated to create that mean 
LF. The associated mean stellar density in this case 
would be given by p, = J < $(M) > dM, or its dis- 
crete summation counterpart (as we shall see in the 
next paragraph, the stability of p* as a function of 
-^simui has been used to define a lower limit for A^simui 
itself). 

At the core of the simulation lies a (pseudo) ran- 
dom number generator. We have tried two different 
generators in order to test the sensitivity of our results 
to the adopted scheme, and found no significant dif- 
ferences in the derived mean overall stellar density, p* , 
as a function of the number of simulations as long as 
the number of simulations is larger than about 1,000 
(see Section pi) . In what follows we have therefore de- 
rived collective values for the LF for 3,000 simulations 
but, evidently, our results are independent of the total 
number of simulations above that minimum number. 
We have found that 3,000 simulations is a good com- 
promise between stability of the simulations, RAM 
memory for array storage, and CPU running time. 
The period (i.e., the number of calls before producing 
correlations) of both random number generators tried 
is, of course, much larger than the number of calls to 
the uniform deviate routine (according to Press et al. 
(1997) the period is on the order of ^ 10^ for our 
adopted generator, see below). For definiteness, we 
have finally adopted the routine "rani" described in 
the last edition of "Numerical Recipes in Fortran", 
from Press et al. (1997). In order to avoid aliasing 
between different "simulations" (as defined in the pre- 
vious paragraph), the "seed" for the random number 
generator is altered between successive simulations. 
Note however that for a given simulation several calls 
to the random number generator are required to pro- 
duce uniform deviates used in the error propagation 
(see Eqn. (Q)). Therefore, the only purpose for up- 
dating the "initial" seed for a given simulation is to 
render successive simulations as differentiated as pos- 
sible. Note also that after providing a seed for "rani" , 
the seed gets altered by "rani" itself, thus taking best 
advantage of the large period in the generator. 

In order to simulate true observational errors, we 
have adopted Gaussian deviates derived using the 
routine "gasdev", also by Press et al. The "seed" 
for the first call to "gasdev" is generated congruently 
with the seed for that particular simulation, but al- 
tered internally in "gasdev" for subsequent calls, as 
per the routine described by Press et al. The value 
adopted for the observable is simply given by, e.g., for 



the parallax: 



Pi.J = Pj + f^Pj X Gi 



(4) 



where pj is the (mean) observed parallax for star 
"j" in the sample, with "measurement error" a-p. , Gi 
is the Gaussian deviate (of zero mean and unity vari- 
ance) for simulation "i" , and pi^j is the i-th simulation 
value for the parallax of star j. The same is performed 
for proper-motion, apparent magnitude, and the bolo- 
metric correction (if necessary). The no-errors situ- 
ation is, of course, reproduced when all the G^'s are 
set to zero. 

A subtlety associated to the simulated sample val- 
ues arises because of the general survey restrictions 
TOb < m < TO/ and /^z < p < Pu. In this case, if a 
simulated value falls outside the survey limits because 
of under/over-shoot due to observational errors, that 
object is eliminated from the sample, and its contri- 
bution to the LF is suppressed. 

6. The current-sample WD LF and its uncer- 
tainty 

In this section we present the results from our nu- 
merical simulations, using the simple technique de- 
scribed in the previous section. As explained before, 
a "convergence" criterion for the simulations has been 
used, by adopting the overall luminosity normaliza- 
tion (stellar density) as the prime parameter. Of 
course, other criteria could be used, but the basic 
point here is that the convergence criterion ensures 
that the derived mean LF becomes independent of 
the number of simulations. It is, therefore, the most 
representative value for the LF given a sample and its 
errors. 

Figure || shows the mean LF resultant from 3,000 
simulations, using the errors quoted by LRB88 and 
LDM98 respectively. Since no errors for the proper- 
motion were given, it was assumed a typical value of 
10 mas yr~^ (where 1 mas= 1 milli-arcsec). As we 
shall see, the exact value for the proper-motion error 
is less critical than uncertainties in the other quanti- 
ties, so this is probably a good estimate. The same 
proper-motions were adopted for both studies. A 
few entries missing errors for V magnitude in LDM88 
were given a probable uncertainty of 0.05 mag, and 
the same magnitude errors were assumed for LRB98, 
which is perhaps an overestimation of LRB98's pho- 
tometric errors, but it is not inconsistent with their 
statement that their "photometric uncertainties are 



3%". The main differences between tliese two stud- 
ies comes from different parallaxes and their un- 
certainties, due to improved parallax series using 
more plates, and the use of CCDs for some of them, 
slightly different optical photometry, also improved 
by the use of digital detectors, and different bolo- 
metric corrections and estimated uncertainties com- 
ing from improved stellar interior and atmospheric 
models and the addition of near-IR photometry used 
in the LRB98 study. We have adopted a flat uncer- 
tainty of 0.14 mag in the BC for the LDM88 sam- 
ple; This is the value they quote for the difference 
between two possible model BC corrections. This is 
probably an overestimation of their true BC errors, 
which applies only to the region of overlap where the 
comparison between different models was done, but 
it provides an upper boundary to their BC errors. 
For the LRB98 sample, we adopt their quoted BC 
uncertainties based mostly on uncertainties in fitting 
their model atmospheres to the broad-band optical 
and near-IR colors. 

Figure y shows a comparison between the no-errors 
WD LF (from Fig. ([|)) and the LF derived using our 
full error simulation (from Fig. (0)). It is apparent 
that, while the mcan-over-simulation WD LF is not so 
different (at least for the brighter bins) in both cases, 
the error bars are quite a bit larger at the fainter bin 
in the Monte-Carlo simulations. As a result, LRB98's 
quoted disc age of 8 Gyr has to be taken with caution 
and, in fact, their data does not rule out the possibil- 
ity of an older disk, with an age as large as 10 Gyr. 
It is also apparent from Fig. that the simulations 
indicate a rather long tail to fainter luminosities in 
comparison with the abrupt decrease seen in the no- 
error calculation. This has, of course, important im- 
plications for the interpretation of the LF in terms of 
a well defined finite age for the disk. Actually, as we 
shall see, this long tail is produced by a few large- 
error excursions in the simulation leading to a biased 
mean LF. Another important point is the exact break 
(if any!) in the LF. This has been an outstanding is- 
sue over the years (c.f. the extensive discussion of this 
on LDM88 or LRB98), the main difficulty here com- 
ing from the fact that the position of the luminosity 
break depends on the exact positioning of the bolo- 
metric magnitude bin, which in the classical 1/I4nax 
method is fixed arbitrarily. Both of these issues have 
motivated us to further explore the behavior of the 
WD LF, specially at the faint end. 

We have devised a simple new version of the 1/ Knax 



method which renders the resulting LF independent 
of the bin positioning. For this purpose, we have 
adopted the same computational scheme of the tra- 
ditional Monte-Carlo method described above, but in 
such a way that we have a "moving box" over abso- 
lute magnitude. The main difference here is that in 
the classical method, the position of the absolute (or 
bolometric) magnitude bins are pre-specified a pri- 
ori, whereas in this new moving box method, we only 
specify the bright and faint absolute (or bolometric) 
magnitude limits, and an arbitrarily large number of 
steps between them (so that the resulting LF appears 
as a continuous function, rather than discrete as in the 
classical case -an important feature if one is trying to 
look for structure in the LF). Because we still need 
the absolute normalization of the LF, the "box" has 
to be integrated over a pre-specified bin width, but 
the position of the luminosity bins themselves can be 
defined over an arbitrarily fine grid. Because the step 
between successive boxes might be smaller than the 
bin width, errors from bin to bin are not totally in- 
dependent, and are thus highly correlated. For this 
reason computing an error at each box position on 
the absolute magnitude grid is less meaningful than 
for the case of the traditional scheme (we can how- 
ever, compute other statistical properties of the LF, 
see below). 

Figure |j shows a comparison of the Monte-Carlo 
LF computed on discrete intervals, and using our 
moving box, both for the LRB98 data set. If, in- 
stead of using the mean LF, we adopt the median 
over the simulations, we reproduce a sharp decline in 
the LF at the faintest bins. Indeed, the quartiles, also 
shown on the figure, indicate a rather tight distribu- 
tion in comparison with the standard deviation from 
the 3,000 simulations. This indicates, in turn, that 
the distribution of predicted LF values, at a given lu- 
minosity, might be quite skewed. If this is the case, 
then, extreme care has to be taken when computing a 
"mean" LF considering the observational errors: One 
must choose an indicator that resembles that of the 
most representative LF value for a given observational 
data set. From this figure we can also see that the me- 
dian LF compares very well with the LF derived using 
no observational errors at all, except in that the latter 
implies a fall-off at slightly brighter luminosities. The 
extreme skewness of the LF is clearly demonstrated 
in Fig. ra where the histogram of LF values is shown 
as a function of bolometric magnitude. We see that, 
as we go to fainter and fainter bins (solid points). 



the LF becomes more and more strongly peaked at 
lower density values, with outliers to large stellar den- 
sity. This is mostly due to large-excursion outliers 
which turn objects intrinsically dimmer in the sim- 
ulations, and thus closer. As a result, the volume 
sampled decreases, and the stellar density increases. 
Another way of appreciating this effect, is to look at 
the stellar density as a function of simulation. Fig- 
ure P shows the predicted overall stellar density as a 
function of simulation, for the 3,000 simulations de- 
scribed above, and for the LRB98 dataset. We see 
that, in comparison with the no-errors predicted den- 
sity, there is a small, but appreciable, over-density 
in the Monte Carlo simulations, due to the sampled 
volume effect just mentioned. As pointed-out before, 
the direct density leads to p = 2.54 x 10~^ stars/pc^, 
whereas a simple fit to the data on Fig. g indicates a 
stellar density of p = (2.59 ± 0.20) x 10"^ stars/pc^ 
(standard error). 

We should emphasize that our derived WD LF is 
still dependent upon bin size. Indeed, the sample 
is still too small to be used as an effective indica- 
tor of different star formation episodes in the disk. 
This is clearly demonstrated in Fig. 0, where we have 
produced three LFs by changing the bin size from 
0.75 mag to 0.25 mag in steps of 0.25 mag, and using 
a moving box with size half that of the bin size, i.e., 
of the same order of the Nyquist frequency for the 
chosen bin size (this was done to avoid strongly cor- 
related errors from bin to bin). As it can be seen from 
this figure, the exact position of the sharp fall-off in 
the LF depends slightly on the bin width chosen, and, 
furthermore, the "wiggly" features in the LF appear 
only for the smallest bin size, indicative of the onset of 
large statistical fluctuations due to the small samples 
concerned (even at relatively bright magnitudes). 

7. The effect of observational errors on the 
WD LF 

In this section we empirically explore the effects of 
different observational error on the derived WD LF. 
We start always from the same sample (for definite- 
ness the LRB98 dataset), but we fudge their errors to 
different amounts in order to understand the behav- 
ior of the 1/Vmax method, and the importance of the 
different sources of errors, on the resultant LF. 

We begin with a WD LF derived using the moving 
box approach described in the previous section. This 
would be the "true" WD LF for this data set (com- 



puted on a continuous set of bins), if there were no 
observational errors. We then add, separately, errors 
in bolometric corrections, magnitudes, parallax, and 
proper-motion, and compare these LFs with that de- 
rived assuming all errors are zero. The outcome of 
these simulations is a prescription, for observers (and 
theorists as well!, see below), as to what parameters 
are more critical, and should thus be refined further. 

Figure plots the quoted observational errors from 
the LRB98 data as a function of apparent V mag- 
nitude. This plot gives us an idea of the range of 
uncertainty in the observational quantities. It is ap- 
parent that the uncertainty in bolometric corrections 
(middle panel, < ctbc >= 0.08 ± 0.02 mag) is a 
lot larger than that of the direct photometry (up- 
per panel, < ay >= 0.028 ± 0.004 mag). This is 
an important point because the WD LF is equally 
sensitive to uncertainties in the BC and apparent 
magnitudes. As for the parallax errors, we have 
< cr^r >= 4.3 ± 0.7 mas. 

Figure ^ shows, as a solid line, the (continuous) 
mean WD LF derived from the LRB98 data by as- 
suming that the uncertainties in all observables are 
zero (this LF "looks" different from the one on Fig. ^ 
which also assumes no errors, because the later uses 
the classical l/Knax method, while in the former we 
are using our "moving box" approach which pro- 
duces a continuous LF). We in turn start "adding" 
errors in various parameters, and discussing their ef- 
fect on the derived LF. In the simulations shown on 
the previous sections, we have adopted a flat error for 
the proper-motion of 10 mas yr~^. The exact value 
adopted is not critical to the resultant LF. Indeed, 
a value 3 times as as large as the assumed one does 
not produce any signiflcant differences in the derive 
LF. Only by the time the errors have gone up to as 
much as 100 mas yr~^ the LF starts showing the ef- 
fects of these errors. Furthermore, the effect only 
appears as a scale (or normalization) factor in the 
overall LF (see upper panel on Fig. ||), and not as 
a signiflcant change of shape on the LF, in marked 
contrast with the effect of errors on apparent magni- 
tude and bolometric corrections (see discussion in the 
next two paragraphs, and the middle and lower pan- 
els of Fig. ^. According to Dawson's analysis of the 
LHS catalogue (Dawson 1986), an intercomparison of 
Luyten's proper-motions and those derived from the 
much more accurate USNO parallax program indi- 
cates that the rms error for a single star in the LHS 
catalogue is 29 mas yr"-'^, in agreement with Luyten's 



own estimates. From about 50% of the sample pub- 
lished so far by Wroblewski & Torres (1990, and sub- 
sequent papers), Costa (2000) finds a sHghtly larger 
dispersion of 38 mas yr~^ (1,262 stars) between their 
proper-motions and those from Luyten for his LTT 
sample (although there is a contribution to this dis- 
persion from their own measurement errors which is 
in the range 5-25 mas yr~^), while Ruiz et al. (2000) 
obtains a dispersion of 36 mas yr~^ (23 stars, with 
similar internal errors as those of the Costa's sample) 
as derived from their respective proper-motion sur- 
veys in selected areas of the southern sky, which has 
recovered many (and added a few new) of the large 
proper-motions stars from the LHS catalogue. We 
thus conclude from these simulations that the current 
proper-motion uncertainties do not contribute signif- 
icantly to the overall LF error. 

Bolometric corrections are, of course, a critical step 
when deriving a WD LF that can be compared with 
theoretical models. Unfortunately, errors in this pa- 
rameter are still rather large (see middle panel on 
Fig. ^, and do produce a large negative impact on 
the accuracy of the derived LF. The middle panel on 
Fig. g shows the effects on the WD LF of uncertain- 
ties amounting to 0.05 mag (red line) and 0.10 mag 
(green line) on the BC. 

Uncertainties in apparent magnitude have the same 
effect on the LF as do uncertainties in the BC, but 
their observational errors are a lot smaller, and there- 
fore do not play an important role in the final WD LF 
(see upper panel on Fig. ||). Parallax uncertainties 
also turn out to be relevant. The bottom panel on 
Fig. P shows the effects of 1 mas (red line), and 3 mas 
(green line) parallax errors on the derived LF. From 
these plots we can clearly see that, for a "typical" un- 
certainty of ctbc ~ 0.1 mag and cttt = 4 mas on these 
two quantities, the contribution to the "smearing" of 
the WD LF is, at present, equally represented by er- 
rors on the trigonometric parallaxes and bolometric 
corrections (all the blue lines on Fig. (g|) indicate the 
contribution of the true errors from proper-motions, 
BC and parallaxes, respectively, to the resultant LF). 

The explanation for the markedly different behav- 
ior of the derived WD LF upon errors on proper- 
motion vs. errors in bolometric corrections and par- 
allaxes is actually easy to understand: As it can been 
seen from Eqns. (1) through (3), the contribution to 
the luminosity function depends on proper-motions, 
parallaxes, and apparent magnitudes, whereas the 
placement in luminosity of a given observed data 



point depends exclusively on the bolometric correc- 
tion, apparent magnitude and (the logarithm of) the 
parallax, but is independent of the object's proper- 
motion. Therefore, as stated before, errors in proper- 
motion will only impact upon the scale of the LF, and 
will not displace points in luminosity, while uncertain- 
ties in the other observables would impact both, the 
LF normalization, and the actual luminosity where 
that object is contributing to the overall LF. 

In terms of the extent to which the sample analyzed 
is complete, our simulations indicate that, with the 
adopted survey boundary constraints, it is probably 
incomplete. Figure (Eol) (upper panel) shows the clas- 
sical estimator < y/Vinax > (which should approach 
the value 0.5 for a complete sample) as a function of 
the simulation ID for the LRB98 sample. As it can 
be seen, no simulation brings the above value closer 
to 0.5. Indeed, for large-excursion errors in which 
one or two objects fall outside the survey boundary, 
the value of < V/VmnK > decreases, leading to two 
parallel sequences to the main set of values. The 
mean error on < y/Vmax > remains, however, con- 
stant and quite small (lower panel on Fig. (|lC|)). We 
can easily explore the effect on the derived LF and 
sample incompleteness due to potentially erroneous 
survey boundaries with the aid of our simulations. In 
Fig. ( pT| ) we show the effect of introducing variations 
in the survey limiting magnitudes, and in the proper- 
motion limits. As mentioned in Section we have 
adopted VJ, — +1, Vj = -1-19, fj,i = 0.8 arcsec yr~-^, 
Hu = 10.0 arcsec yr~^. We find that the derived LF is 
insensitive to the value of Vf,, and that even by adopt- 
ing the extreme case Vj, — +5 the LF does not change 
at all. Also, surprisingly, by adopting a very conser- 
vative cut Vf = -1-18, the LF is only altered mildly. 
Proper-motions do have, however, an important role 
in the sample selection, and in the resultant LF (as 
found also by the simulations from Wood and Oswalt 
(1998) and Garcia-Berro et al. (1999)). While a cut 
at ^u — 5.0 arcsec yr~^ does not change apprecia- 
bly the LF, a lower value oi ^u = 2.0 arcsec yr~^ 
does change the shape of the LF appreciably. The 
only remaining source of uncertainty, i.e., the lower 
proper- motion limit, does also seem to have a big in- 
fluence on the derived LF but only for Mtoi < -1-15.0 
(see Fig. (pT|)), and only when adopting a very con- 
servative m = 1.0 arcsec yr~^. Figure ( pl| ) also 
shows the encouraging news that the exact break at 
the faint end of the WD LF is not extremely sensi- 
tive to the survey boundary and/or incompleteness 



effects. From a comparison of Luyten's catalogued 
stars and their newly discovered large proper-motion 
stars in selected areas of the Southern sky, Ruiz et 
al. (2000) find that the LHS catalogue might be 
actually incomplete in a more severe way than pre- 
viously thought, at least in the Southern sky. For 
example, while they corroborate that the LHS cat- 
alogue is incomplete for fi < 0.8 arcsec yr~^, they 
also find that incompleteness sets for ttir > 14, i.e., 
several magnitudes brighter than claimed by Luyten, 
although the overall number of stars involved in the 
comparison is small (~ 50). This might not neces- 
sarily apply directly to the Northern sky sample ana- 
lyzed here, since it is known that the Southern plates 
used by Luyten (those of the Bruce proper-motion 
survey) were shallower {rripg^.^^ ^ 15.5 — 16.0) than 
its Northern counterpart (from the Palomar proper- 
motion survey, with rUpg^.^ ~ 21.2). One must note 
that the completeness test for the noise-free classi- 
cal case also shows signs of incompleteness, having a 
value of < y/Knax >= 0.367 ± 0.046 for the LRB98 
dataset (hence the suggestion of a significant incom- 
pleteness or erroneous survey boundaries is nothing 
particular to our simulations!). 

From the preceding simulations, it is interesting to 
notice that the only way we can increase the value 
of < y/Kiiax > is by adopting a smaller //„. A 
value of ^„ = 2.0 arcsec yr~^ imphes < ^/Vmax >= 
0.447 ± 0.051, whereas /i„ = 1.5 arcsec yr^^ leads 

< WKnax >= 0.594 ± 0.062 (see Fig. (^). By com- 
parison, the adopted standard survey boundaries have 

< V'/Kiax >= 0.366 ± 0.046. If the incompleteness 
is due to a poor definition of the survey boundaries, 
then this suggests that Luyten might have actually 
missed some of the large (/i > 1.5 — 2.0 arcsec yr~^) 
proper- motion stars, a fact also stated by Dawson 
(1986) in connection with the luminosity function 
of halo stars, as derived from the LHS catalogue. 
Our simulations also imply that he only needed to 
have missed stars moving with proper-motions larger 
than 1.5 to 2.0 arcsec yr~^, but smaller than 4.0 
to 5.0 arcsec yr~^, since the derived LF (and the 
corresponding completeness test) is insensitive to fi^ 
for ^u > 5.0 arcsec yr"-'^. Even though the LF 
is influenced heavily for Mt,oi < -1-15.0 by adopting 
^i = 1.0 arcsec yr~^, < V/V^ax > does not change 
substantially in this case (0.372 ± 0.050), indicating 
that the source of incompleteness is due to, both, faint 
and large proper-motion stars missing from the cur- 
rent sample. 



8. Discussion and Conclusions 

As was explained towards the end of Sect, g, sim- 
ulated sample points whose observables fall outside 
the survey limits are eliminated from the overall sam- 
ple. This procedure might be called into question 
because the resulting WD LF has a different num- 
ber of stars than the number contained in either the 
LDM88 or LRB98 samples. Indeed, this strategy can 
only reduce the number of objects in the realization, 
whereas in the real world, objects can be added to 
the sample as well, since the same observational er- 
rors will occasionally add stars to the sample that 
were originally outside the survey limits. The non- 
conservation of sample data points is clearly seen in 
the second, parallel sequence in < y/Vmax > that 
appears in Figs.Qr 



and 12 



The question is, then, to 
what extent the simulations displayed in Fig. 2, which 
show a long tail of intrinsically faint stars with small 
implied observational errors, might be an artifact of 
the non-conservation of the total numbers of stars in 
the (simulated) sample? 

In Sect, n we have interpreted the small value for 
< V/Vmax > (see also Fig.|lO|) as a reflection of an 
incorrect deflnition of the survey limits, and we point 
out that the second sequence in < ^/Vmax > is due 
to large-excursion errors in which one or two objects 
fall outside the adopted survey boundary. This is the 
same explanation we advanced for the faint tail of the 
WD LF in Sect. g. But, one might also question this 
result, and wonder whether the fact that this quan- 
tity departs from the expected value of 0.5 may not be 
a reflection of an incorrect treatment of these large- 
excursion errors, as suggested in the previous para- 
graph, rather than reflecting an incorrect treatment 
of the survey limits. 

To elucidate these important questions we run a 
few more simulations where we discarded the entire 
realization of the WD LF if the number of objects 
was not conserved, and proceeded to the next real- 
ization. In this case the simulated samples were, by 
construction, always similar to the input sample, save 
for the exact values of the observables which departed 
form the input values by an amount specified only by 
the adopted errors, but still within the adopted sur- 
vey boundaries. The results of these simulations are 
shown in Fig. n3 (upper panel), which compares the 
LF derived in the case where the sample is not con- 
served (i.e., the same LF as for the LRB98 dataset 
shown in Fig. 0, derived from 3,000 simulations) with 



10 



the LF obtained in the case when the sample is con- 
served (2,557 simulations out of 3,000 initial simula- 
tions, i.e., some 15% of the simulations lost to ex- 
tremely high-residual excursions). As it is obvious 
from this figure, the resultant LFs are quite similar, 
and both exhibit the same tail to faint luminosities. 
This implies that the large-excursion simulated ob- 
servables responsible for the faint tail of the LF are 
not large enough to alter the input sample signifi- 
cantly (therefore avoiding the problems mentioned in 
the first paragraph above) , while having an important 
effect on the faint portion of the LF itself due to the 
highly skewed distribution of LF values at a given lu- 
minosity shown in Fig. |5|. The overall completeness 
factors do not change either, having in both cases the 
value < T^/Vmax >= 0.366 ± 0.046 (and quite close to 
the no-error case with < 1//Vmax >= 0.367 ± 0.046), 
thus alleviating the concerns expressed in the previ- 
ous paragraph. 

In light of the above results, it is interesting to ex- 
plore in some detail the selection effects acting upon 
the derived simulated LFs and inclusion/exclusion of 
objects near the survey limits. For this purpose, we 
calculated the number of times that the selection ef- 
fect was either proper-motion or apparent magnitude, 
according to the first and second terms respectively 
of Eqs. H and g, as a function of bolometric magni- 
tude. For the bolometric WD LF exhibited in the 
upper panel of Fig. O, r^in was always determined 
by the proper-motion, while the value of the Tmax 
was determined by a mixture of both proper-motion 
and apparent magnitude. For A^obj objects in the 
input sample, and Nsimui simulations, the maximum 
number of times that a criteria could be used is just 
A'obj X -^simui- In the lower panel of Fig. Owe show 
the histogram of selection criteria for r^nax for 43 ob- 
jects (LRB98 dataset), and 2,557 effective simulations 
(keeping the number of objects fixed, see previous 
paragraph). As it can be seen from the figure, in the 
vast majority of the simulations, the primary selection 
criteria in the whole range 13 < Afboi < 17.5 is deter- 
mined by the object's proper-motion, and not by its 
apparent magnitude, a result already found by Wood 
and Oswalt (1998) and Garcia-Berro et al. (1999, 
see related discussion below). Additionally, it was 
demonstrated in Sect. and Fig. g (top panel), that 
errors in proper-motion will only impact upon the 
scale of the LF, but will not displace points in lumi- 
nosity, while uncertainties in the other observables (in 
particular apparent magnitude and bolometric correc- 



tions) will impact both, the LF normalization, and 
the actual luminosity where that object is contribut- 
ing to the overall LF. We therefore can reasonably 
suggest that the true significance of the by and large 
proper-motion selection criteria at faint magnitudes, 
intermingled with a much less often magnitude selec- 
tion criteria, is responsible for keeping the number of 
objects fixed within the survey boundaries (proper- 
motion selection, changes LF normalization but not 
the object's luminosity) while producing a long tail to 
faint magnitudes (magnitude selection, changes nor- 
malization and luminosity, but only in a small fraction 
of the simulations, see Fig. p^ ). 

We can therefore conclude that our simulations in- 
dicate that LRB98's data, when properly accounting 
for observational errors, does not rule out a disk with 
an age as large as 10 Gyr. This is good news because 
previous studies that find ages of 8 Gyr or younger 
using similar datasets are difficult to reconcile with 
an halo age of 15 Gyr (inferred from old globular 
clusters) given that Galactic formation and chemical 
evolution models suggest a delay of, at most, 3 Gyr 
between the onset of star formation in the halo and 
in the local disk (Wood and Oswalt 1998). On the 
other hand, a value as large as 13 Gyr, found by 
Garcia-Berro et al. (1999), can probably be ruled 
out. Garcia-Berro et al. attribute this large derived 
age to the effects of the scale height, but we should 
note that our study is restricted to the faintest por- 
tion of the LF (with rmax from Eq. (|l|) always smaller 
than 50 pc, or about 1/3'''* of a scale height), where 
the effect of the WD scale-height is irrelevant because 
of the local nature of the sample. We note that Galac- 
tic open clusters can also set constraints on the age of 
the disk. NGC 6791 is the oldest known open cluster 
with a well-determined age, in the range 7 to 10 Gyr 
(Tripicco et al. 1995), albeit Scott el al. (1995) find 
that the kinematics (space velocity) for this object is 
somewhat peculiar. The extreme 12 Gyr age of Berke- 
ley 17, believed to be one of the oldest disk clusters 
(Phelps 1997), has been recently revised by Carraro 
et al. (1999) using near-IR photometry, leading to a 
younger age of 8-9 Gyr. Also, Jimenez et al. (1998), 
using Hipparcos data, have found an upper limit for 
the age for the disk field population in the solar neigh- 
borhood of 11 ± 1 Gyr (see also Bertelli et al. 1999), 
which would be in agreement with our revised (older) 
age from the WD LF. Also, we find that current obser- 
vational uncertainties and sample sizes do not allow 
us to establish the existence of small scale features 
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in the WD LF which could be indicative of different 
episodes of star formation in the disk. This could only 
be alleviated by dramatically increasing the currently 
small samples, as also emphasized by the simulations 
performed by Wood and Oswalt (1998), and Garci'a- 
Berro et al. (1999). 

Both Wood and Oswalt (1998) and Garcia-Berro 
et al. (1999) have found that the primary selec- 
tion criteria at low luminosities is the proper-motion. 
In our simulations we see a related effect, where a 
larger proper-motion uncertainty (which affects the 
selection criteria) leads to a change of the normal- 
ization of the LF, whereas its overall shape does not 
change dramatically (see Fig pi, upper panel) . By us- 
ing the < V^/Knax > completeness criteria, we also 
find that the LRB98 sample seems to be missing faint 
{Mboi > -1-15.0), large proper-motion (/x > 2.0 arc- 
sec yr~^) stars, and that the sample is only com- 
plete for II < 1.5 — 2.0 arcsec yr^^. However, we 
find that the precise luminosity break at the faint 
end of the WD LF is not extremely sensitive to the 
survey boundary and/or incompleteness effects (see 
Figure (0)). 

In summary, we have found that most of the 
current uncertainties in the observational WD LF 
come from uncertainties in bolometric corrections and 
in parallaxes, while photometry and proper-motions 
play a minor role. Although this effect might be cap- 
tured in Fig. @, which simply displays the distribution 
of errors in the bolometric corrections and parallaxes 
(and which, of course, does not require any of the sta- 
tistical discussion in the rest of the paper) , the impact 
of these uncertainties upon the Monte-Carlo derived 
WD LF for 3,000 simulations is fully shown in Fig. O. 
This last figure clearly shows that, refinements on the- 
oretical models (such that ctbc < 0.05 mag) and par- 
allaxes (with (Ttt < 1 mas), as well as larger samples 
(A^samp ~ 200, see Wood and Oswah 1998), should be 
primary goals in order to produce a better luminosity 
function for white dwarfs. 

Our Monte-Carlo simulations using the l/I^nax 
method can be recreated for any other proper-motion 
and/or magnitude selected samples - a simple ASCII 
input table with mean values and probable errors is 
all what it is needed. All the software and help on 
how to use it is available from RAM. 
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Fig. 1. — A comparison of the LDM88 (solid cir- 
cles) and LRB98 (solid squares) WD LF. A bin width 
of 0.5 mag, and the same bin centering adopted by 
LDM88 and LRB98 have been used. A shift of 
+0.1 mag in the magnitudes has been applied to the 
LRB98 data in order to avoid crowding. The up- 
per panel is for A/y, while the lower panel is for 
M})oi , adopting the bolomctric corrections applied by 
LDM88 and LRB98 respectively. Errors bars indicate 
Poisson error exclusively. This figure shows, basically, 
our advance in 10 years of optical & IR photometry, 
parallaxes, bolometric corrections and interior physics 
and stellar atmosphere modeling for these stars, as re- 
flected upon the LF. 

Fig. 2. — Similar to Figure nl but for the mean- 
over-simulations WD LF for 3,000 Monte-Carlo sim- 
ulations adopting the errors quoted by LDM88 and 
LRB98. By comparing with Fig. nl it is apparent 
that the Monte-Carlo errors are smaller than the Pois- 
son errors at brighter magnitudes where the sample is 
larger and less subject to observational errors. How- 
ever, we can also see that Poisson errors underesti- 
mate the true expected LF uncertainties at the fainter 
bins. The mean value for the LF are, however, quite 
similar in both cases. 

Fig. 3. — Bolometric WD LF. Similar symbols as for 
Fig. 0. In this case, LRB98 has been shifted by +0M 
in \og{L / Lq) in order to avoid crowding. The upper 
panel is for no-observational errors, while the lower 
panel shows the results for the mean LF after 3,000 
simulations. The lines indicates the latest theoretical 
WD LF published by Benvenuto and Althaus (1999), 
based on carbon-oxygen core WDs, including the re- 
lease of latent heat during crystallization (Salpeter 
IMF, constant star- formation rate over the age of the 
disk). The dotted line is for a 6 Gyr disk, while the 
dot-dashed is for a 10 Gyr disk. Theoretical LFs have 
been normalized to a total p — 3.39 x 10"'^ stars/pc^. 
The upper panel shows that, as found by LRB98, the 
WD LF is not inconsistent with a disk age of about 
8 Gyr. However, in the lower panel, our Monte-Carlo 
approach indicates that the large error-bars at the 
faint end could allow for an older disk given current 
observational uncertainties and sample sizes. 

Fig. 4.— Bolometric LF for the LRB98 dataset. The 
solid squares with error bars indicate the LF using 
the Monte-Carlo mean LF on discrete 0.5 mag bin 
intervals (same as in Fig. 0, lower panel), while the 



open squares reproduces the LF derived in the case of 
no errors, shifted by -1-0.04 (same as the solid squares 
in Fig. m . The solid line shows the mean LF using our 
moving-box approach. In this case, the LF has also 
been integrated over a 0.5 mag bin. The dashed line 
indicate the median over simulations LF from the very 
same simulation that generated the plotted mean LF, 
while the dotted lines indicate the lower 25% quartile 
and the upper 75% quartile on the distribution of LF 
values as a function of luminosity. The big difference 
between the mean and median LF at faint magnitudes 
indicates a highly skewed distribution of LF values, as 
it is indeed found (see Fig. g). Surprisingly, we can 
also see that the median LF approaches very well the 
LF derived in the case of no errors. 

Fig. 5. — Distribution of Monte Carlo LF values as 
a function of bolometric luminosity for the LRB98 
data set in 0.5 mag bins. The plot shows the evo- 
lution of the histogram of LF values as a function 
of bolometric magnitude. The LF has been derived 
from 3,000 simulations leading to the mean and me- 
dian LF for LRB98 shown in the previous figures. 
The symbols are, in decreasing luminosity, as fol- 
lows: Open circles for A'/boi = +14.5, open squares for 
-Mboi — +15.0, open triangles for Mboi = +15.5, open 
stars for il/boi — +16.0, filled circle for il/boi = +16.5, 
filled squares for A/boi = +17.0, and filled triangles 
for Mboi = +17.5. 

Fig. 6. — Predicted WD space density as a function 
of simulation, for the 3,000 simulations leading to the 
LF on the previous figures (dots). The horizontal 
solid line indicates the density derived in the classi- 
cal no-errors method. Changes in the observables al- 
lowed within their uncertainties lead to a "spill-over" 
to adjacent luminosity bins that increase the effective 
density by a small fraction due to a volume sampling 
effect (see text). Although there is a systematic shift, 
the dispersion across simulations is large enough that 
the Monte-Carlo mean density and the direct density 
are almost the same within the dispersion of the for- 



Fig. 7. — Continuous median WD LF from the 
LRB98 data set adopting three different bin widths of 
0.75 mag (dashed line), 0.5 mag (solid fine, adopted 
throughout this paper), and 0.25 mag (dot-dashed 
line). As it can be seen, rapid changes in the LF ap- 
pear only for the smallest bin width due to the small 
number of objects per bin. At the faint end, the small 
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bin- width effect is seen as an uncertainty in the exact 
fall-off luminosity. 

Fig. 8. — Quoted observational errors (dots) as a 
function of apparent V magnitude for the LRB98 sam- 
ple. Proper-motion errors have been assumed to be 
10 mas (see text), and are thus not plotted. The up- 
per, middle and lower panels are for uncertainties in 
apparent magnitude, BCs, and trigonometric paral- 
laxes respectively. The horizontal line indicates the 
straight mean error for all quantities. 

Fig. 9. — Linear continuous (using our moving box 
approach) WD LF as a function of uncertainties in 
proper- motion (upper panel), BC (middle panel), and 
parallax (lower panel). In all panels, the black solid 
line is for no errors while the blue solid line is for the 
quoted (true) errors. In the upper panel the red line is 
for a proper-motion uncertainty of 30 mas yr~^ (more 
representative of the LHS catalogue) , while the green 
line is for an extreme error of 100 mas yr^^, larger by 
a factor of three than the expected errors in the LHS 
catalogue. As it can be seen from the plot, proper- 
motions are not a significant source of uncertainty on 
the derived LF — they only affect the bin-to-bin nor- 
malization of the LF, but do not broaden the lumi- 
nosity distribution. In the middle panel, the red line 
is for an uncertainty of ctbc = 0.05 mag, while the 
green line is for an uncertainty of ctbc — 0.10 mag. 
In the lower panel, the red line is for a parallax un- 
certainty of only 1 mas, while the green line is for an 
uncertainty of 3 mas. The largest source of uncer- 
tainty in the present WD LF is found to come from 
uncertainties in both, the bolometric corrections and 
the trigonometric parallaxes. 

Fig. 10. — Completeness fraction (top panel), mea- 
sured in terms of < V/X^nax > as a function of simu- 
lation, for the 3,000 simulations resulting in the bolo- 
metric LF of Fig. g. The lower panel shows the ex- 
pected mean error in < y/Vmax > for the respective 
simulations. The horizontal line indicates the values 
for the classical error-free case with < F/Vmax >= 
0.367 ± 0.046. A truly complete sample should have 
< WKnax >-^ 0.5 (Schmidt 1968,1975, Felten 1976) 
whereas the LRB98 dataset has < l^/l^nax >< 0.38, 
indicating that either the sample is somewhat incom- 
plete, or the survey boundaries are erroneous. 

Fig. 11. — Linear continuous median WD LF for dif- 
ferent survey boundaries. The solid line is for the 



standard survey boundaries (dashed line on Fig. ^, 
the dashed line for /i; — 1.0 arcsec yr^^, the dot- 
dashed line for yu„ — 2.0 arcsec yr~^, and the dotted 
line is for the extreme case /i„ = 1.5 arcsec yr~^. As 
can be shown, the lower proper-motion limit does not 
have any effect on the fall-off luminosity, but it has 
a great influence on the LF for Mi,oi < -f 15.0. The 
large proper-motion cut, however, influences both the 
bright & faint portions of the LF and, in addition, 
changes the completeness of the sample (see Fig. O.) 

Fig. 12. — Similar to Fig. dlQ), but with different val- 
ues for the survey boundaries (see text). Green dots 
are for /i„ = 2.0 arcsec yr~^ (with < V/Vma^ >= 
0.447 ± 0.051), while red dots are for /!„ = 1.5 arc- 
sec yr^^ (with < 1//T4iax >= 0.594±0.062). This sug- 
gests that the data is complete ( < F/Knax >= 0.5) 
only for /i„ > 1.7 arcsec yr~^. 

Fig. 13. — Top panel: Mean discrete bolometric 
WD LF for the LRB98 dataset for 3,000 Monte-Carlo 
simulations discarding individual simulated objects as 
their errors move them out of the survey boundaries 
(solid squares - same LF as solid squares in lower 
panel of Fig. g), or by discarding the simulated LF 
altogether from the computation of the mean LF if 
the simulated value falls off the adopted boundaries 
(open squares, 2,557 effective simulations). As it can 
be readily seen, the long tail to faint luminosities re- 
mains intact even if we exclude the influence of ex- 
treme simulated outliers while keeping the sample size 
fixed. Lower panel: Log of the number of times that 
either a proper-motion (solid line) or apparent mag- 
nitude (dashed line) criteria is adopted to calculate 
''max from Eq. [^ to compute the (open squares) LF 
displayed in the top panel. As it is obvious, the pri- 
mary selection effect at faint magnitudes is proper- 
motion (see text for significance of this). 
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